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Abstract 

ON 

. We combine the analyses for flavor changing neutral current processes and dark matter solutions 
in minimal-type supersymmetric grand unified theory (GUT) models, SO(10) and SU(5), with 
a large B s -B s mixing phase and large tan/3. For large tan/3, the double penguin diagram 
dominates the SUSY contribution to the B s -B s mixing amplitude. Also, the By(B s — > /i/i) 
constraint becomes important as it grows as tan 6 (3, although it can still be suppressed by large 
pseudoscalar Higgs mass m^. We investigate the correlation between B s — > \x\i and the dark 
matter direct detection cross-section through their dependence on m^. In the minimal- type of 
SU(5) with type I seesaw, the large mixing in neutrino Dirac couplings results in large lepton 
flavor violating decay process r — > /ry, which in turn sets upper bound on tua- In the SO(10) 
case, the large mixing can be chosen to be in the Majorana couplings instead, and the constraint 
from Br(r — ► /ry) can be avoided. The heavy Higgs funnel region turns out to be an interesting 
possibility in both cases and the direct dark matter detection should be possible in the near 
future in these scenarios. 



1 Introduction 



Recently, CDF and D0 collaborations have announced the analysis of the flavor-tagged B s — > 
J/i[>4> decay. The decay width difference and the mixing induced CP violating phase, <f) s , have 
been extracted from their analysis pQ. In the Standard Model (SM), the CP violating phase is 
predicted to be small, 4> s = 2/3 s = 2 arg (—V ts V t * b /V cs V* b ) ~ 0.04. However, the measured values 
of the phase are large: 

s (CDF) G [0.28,1.29] (68%C.L.), (1) 
S (D0) = 0.57t°J°(stat)l°$(syst). (2) 

Combined data analyses including the semileptonic asymmetry in the B s decay indicate that 
the CP violating phase deviates about 3cr from the SM prediction [2J. If this large phase still 
persists in the upcoming results from Fermilab, it implies the existence of new physics (NP) 
beyond the SM and that the new physics requires a flavor violation in b-s transition as well as 
a phase in the transition. 

Supersymmetry (SUSY) is the most attractive candidate to build NP models. As it is 
well known, SUSY models have a natural dark matter candidate which is the lightest SUSY 
particle (LSP). Besides, the gauge hierarchy problem can be solved and a natural aspect of the 
theory can be developed from the weak scale to the ultra high energy scale. In fact, the gauge 
coupling constants of the Standard Model gauge symmetries can unify at a high scale using the 
renormalization group equations (RGEs) involving the particle contents of the minimal SUSY 
standard model (MSSM). This indicates the existence of grand unified theories (GUTs). The 
well motivated SUSY GUTs have always been subjects of intense experimental and theoretical 
investigations. Identifying a GUT model, as currently is, will be a major focus of the upcoming 
experiments. 

The nature of the flavor changing neutral currents (FCNCs) and the CP violating phase 
is very important to test the existence of new physics beyond the standard model. In SUSY 
models, the SUSY breaking mass terms for squarks and sleptons must be introduced, and they 
have sources of FCNCs and CP violation beyond the Kobayashi-Maskawa theory. In general, 
the soft breaking terms generate too large FCNCs, hence flavor universality is often assumed 
in squark and slepton mass matrices to avoid large FCNCs in the meson mixings and lepton 
flavor violations (LFV) [3]. The flavor universality is expected to be realized by the Planck 
scale physics. However, even if universality is realized at a scale such as the GUT scale or the 
Planck scale, non-universality in the SUSY breaking sfermion masses is still generated from 
the evolution of RGEs, and this can lead to a small flavor violating transitions, which could be 
observed in the ongoing experiments. 
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In the MSSM with right-handed neutrinos, the induced FCNCs from the RGE effects are not 
large in the quark sector, while sizable effects can be generated in the lepton sector due to the 
large neutrino mixing angles [I]. Within GUTs, however, loop effects due to the large neutrino 
mixings can induce sizable FCNCs also in the quark sector since the GUT scale particles which 
connect quark and lepton sectors can propagate in the loops [5]. As a result, the patterns of 
the induced FCNCs highly depend on the unification scenario and the heavy particle contents. 
Therefore, it is important to investigate FCNC effects to obtain a footprint of the GUT models. 
If quark-lepton unification is manifested in a GUT model, the flavor violation in b-s transition 
can be responsible for the large atmospheric neutrino mixing |6j, and thus, the amount of the 
flavor violation in b-s transition (the second and the third generation mixing), which is related 
to the B s -B s mixing and its phase, has to be related to the r — > 717 decay [TJ, El El DI]] for a 
given particle spectrum. The branching ratio of the r — > ^7 decay is being measured at the 
I?-factory, and thus, the future results on LFV and from the ongoing measurement of the phase 
of B s -B s mixing will provide an important information to probe the GUT scale physics. 

In Refs.jHJ HI], two of us have studied the correlation between Br(r —* 717) and <fi s , the 
phase in B s -B s mixing, comparing between SU(5) and SO(10) GUT models, and investigated 
the observational constraints in these models in order to decipher the GUT physics. The flavor 
violation, originating from the loop corrections via heavy particles, can be characterized by the 
CKM (Cabibbo-Kobayashi-Maskawa) quark mixing matrix and the MNSP (Maki-Nakagawa- 
Sakata-Pontecorvo) neutrino mixing matrix, as well as the size of the Yukawa couplings. Since 
the CKM mixings are small, it is expected that the neutrino mixings dominate the source of 
FCNCs at low energy. It is important to know whether the large neutrino mixings originate 
from Dirac-type or Major ana-type neutrino Yukawa couplings. When the large neutrino mixings 
originate from the Dirac neutrino Yukawa couplings in a particular GUT model, the (squared) 
right-handed down- type squark mass matrix, M|, c , as well as the left-handed lepton doublet 
mass matrix, M|, can have flavor non- universality. When the large mixings originate from the 
Majorana Yukawa couplings, the left-handed squark mass matrix, M? , can also have flavor non- 
universality in addition to the other sfermions. In the minimal-type of SU(5) GUT, the large 
neutrino mixings originate from Dirac neutrino couplings provided that there is no fine-tuning 
in the seesaw neutrino matrix. On the other hand, in the minimal- type of SO (10) GUT, the 
large neutrino mixings can originate from Majorana-type couplings. In general, since SU(5) is a 
subgroup of SO(fO), one can construct a model, a non-minimal-type of SU(5) GUT, where the 
neutrino mixings originate from Majorana-type couplings. Conversely, if we allow fine-tuning 
in the Yukawa coupling matrices, Dirac neutrino Yukawa couplings can be the source of the 
large mixings even in the SO(fO) model. Actually, there is a little ambiguity to determine 
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the minimal SU(5) or SO(10) GUT model because the very minimal versions of the GUT 
models have problems with phenomenology and a slight modification is needed. (That is why 
we call our models "minimal-type".) Here, we call a typical boundary condition as minimal- 
type SU(5) GUT condition when the off-diagonal elements of M? c and M| are correlated due 
to the Dirac neutrino couplings in the model. Another type of boundary condition where 
the Mq is correlated to Mj~ c ^ c and M| due to the Majorana coupling in SO(10) model is 
called the minimal-type SO(10) GUT boundary condition. The large phase of B s -B s mixing 
in combination with the other flavor violating processes, can tell us which type of boundary 
condition is preferable. 

We analyzed the case of a large tan (3 (which is the ratio of the vacuum expectation values of 
up- and down-type Higgs fields) in the Ref . [IT] . In large tan (3 case, the so-called double penguin 
contribution [T2J [13] can dominate the SUSY contribution to the B s -B s mixing amplitude 
over the box contribution unless the pseudo Higgs field is heavy. When the double penguin 
contribution is enhanced by a smaller pseudo Higgs field mass, the B s — ► /x/x decay [HI |J3] is 
also enhanced close to its experimental bound [J5]. In other words, if the large phase of B s -B s 
mixing originates from the double penguin contribution, the B s —>■ /x/x decay should be observed 
very soon, and therefore it is worthwhile to examine the constraints to see if the large phase 
is really generated from the double penguin contribution. An important constraint to obtain 
a large phase of B s -B s in GUT models comes from the experimental bound of r — > /X7 decay. 
Due to the quark-lepton unification, when the flavor violation of b-s transition is large, the r-/x 
flavor violation is expected to be large as well. This is significant especially in the minimal-type 
of SU(5) model. As a result, a wide region of the parameter space is excluded. This result is 
important to distinguish among the solutions of the dark matter relic density as measured by 
Wilkinson Microwave Anisotropy Probe (WMAP) [16] in the context of SUSY dark matter. 
Since imposing supersymmetric solution to the dark matter content of the universe puts a 
tremendous constraint on the SUSY model parameter space, it is interesting to see whether 
the constraints from the flavor violating processes can be satisfied by the dark matter allowed 
regions. The allowed parameter space, satisfying both of these constraints, can then be probed 
directly at the large hadron collider (LHC). In this paper, we assume the lightest neutralino 
as dark matter in the GUT models and combine the dark matter analyses with the flavor 
constraint analyses. Consequently, some of the solutions for the neutralino dark matter are 
disfavored, and the so-called funnel solution (in which the neutralinos annihilate through the 
heavy Higgs bosons) turns out to be an interesting one for the large phase of B s -B s mixing 
and large tan/3 case considered here. In the solution which satisfies all the constraints, the 
branching ratio of the B s — > /x/x decay is predicted, and we will show that it is in the range 
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to be observed soon. In scenarios with large tan/3, the dark matter direct detection can be 
correlated to the branching ratio of the B s — > fi/j, decay since both are enhanced by a small 
pseudoscalar Higgs mass [17]. We investigate this correlation for the case with large phase of 
B s -B s mixing. 

The paper is organized as follows: In section 2, we describe the FCNC sources in SUSY 
GUT models. The two typical boundary conditions in both SU(5) and SO(10) models are 
considered. In section 3, we describe the SUSY contributions of B s -B s mixing amplitudes, and 
the constraints from the other FCNC modes. We also discuss the solutions of the WMAP dark 
matter relic density which can be allowed in this scheme. In section 4, we show our numerical 
results on both SU(5) and SO(10) GUT models. Section 5 is devoted to conclusion and remarks. 

2 FCNC sources in SUSY GUTs 

In SUSY theories, the SUSY breaking terms can be sources of flavor violations. In general, it 
is easy to include flavor violating terms by hand since the SUSY breaking masses with flavor 
indices are parameters in the MSSM. However, large FCNCs are induced if these parameters are 
completely general [3]. Therefore, as a minimal assumption of the SUSY breaking, universality 
of the scalar masses is often considered. This means that all the SUSY breaking (squared) 
scalar masses are equal to and all the scalar trilinear couplings are proportional to the 
Yukawa couplings (with the coefficients are universal to be A Q ) at a unification scale. 

Even if the universality is assumed, non-universality in the scalar masses at the weak scale 
is generated by the evolution of the theory from the GUT scale down to the lower scale via the 
RGEs. As we have mentioned in the introduction, in the MSSM with right-handed neutrinos 
(iV c ) the induced FCNCs from the RGE effects are not large in the quark sector while sizable 
effects can be generated in the lepton sector due to the large neutrino mixings pi] . The sources of 
FCNCs in this model are the Dirac neutrino couplings. In GUT models, the left-handed lepton 
doublet (L) and the right-handed down- type squarks (_D C ) are unified in 5, and the Dirac 
neutrino couplings can be written as Y U 5N C H^. As a result, non-universality in the SUSY 
breaking mass matrix for D c is generated from the colored-Higgs and right-handed neutrino 
loop diagrams, and flavor violation in the quark sector can then also be generated from the 
Dirac neutrino couplings [5j [6] . 

The light neutrino mass matrix can be written as 

M** = /(A L ) - Y U M R X T (H° U ) 2 , (3) 

where Al is an SU(2)^ triplet, and / is a Majorana coupling ^LLA^. The second term is 
called type I seesaw term [18j. If the type I seesaw term dominates the light neutrino mass, 
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the Dirac neutrino coupling must have large mixings to explain the large neutrino mixings in 
the basis where the charge-lepton Yukawa coupling matrix Y e is diagonal. On the other hand, 
when the first term (triplet term) dominates (type II seesaw |19|), the Majorana coupling must 
have large mixings. Distinguishing these two cases is very important in order to understand 
the source of FCNCs in the GUT models. 

The triplet contribution in the type II seesaw is natural in the framework of SO(10) GUT 
models |20j. In the SO(10) models, all matter species are unified in the spinor representation 16. 
Since the right-handed neutrino is also unified to other matter fields, the neutrino Dirac Yukawa 
coupling does not have large mixings in the minimal-type of SO(10) models, and the proper 
neutrino masses with large mixings can be generated from the Majorana couplings |/LLA^. 
The / coupling is unified to the 16 16 126 term which also includes Dirac Yukawa couplings 
for fermions, and thus the model is predictive [21J. If any of the decomposed fields from 
126 is lighter than the unification scale, the flavor non-universality for squarks and sleptons 
is generated. It is then possible that the loop corrections generate the flavor violation for 
the left-handed quark doublet (Q), the right-handed up- type quark {U c ) and the right-handed 
charged- lepton (E c ), in addition to D c and L. 

We parameterize the non-universality in the squark and slepton mass matrices due to the 
loop corrections as 

M) = mg[l - K F U F dmg(k 1} k 2 , 1)4], (4) 

where F = Q, U c , D c , L, E c . The quantity Kp denotes the amount of the off-diagonal elements 
and it depends on the sfermion species. The unitary matrices Up is equal to the neutrino mixing 
matrix in a limit. We note that the unitary matrices Up should be defined in the basis where 
charged-lepton Yukawa matrix, and down-type quark Yukawa matrix are diagonal in order to 
calculate the flavor violating processes such as r — > /ry, and B s -B s mixing. In the minimal 
SU(5) GUT where only H 5 and H$ couple to the fermions by renormalizable terms, Up>c is 
exactly same as Ul and has large mixing angles, while Uq xjo e° have small mixings relating the 
CKM mixings. In general, fermion mass matrices come from the sum of the Yukawa terms, 
and the equality of Up>c and Ul can be completely broken when there are cancellations among 
the minimal Yukawa term and additional Yukawa terms. Here, we consider a model where the 
(nearly) equality between Up>c and Ul (especially for 23 mixing angle of them) is maintained 
as a "minimal-type" assumption. The assumption is natural if there is a dominant Yukawa 
contribution and corrections to fit realistic masses and mixings are small. In the minimal-type 
of SO (10) model, all Up can have large mixings responsible for the neutrino mixings. The 
detail physical interpretation of this parameterization is given in [TTJ [22]. When the Dirac 
neutrino Yukawa coupling Y v or the Majorana coupling / is hierarchical, we obtain k%, k 2 <C 1 
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and then the 23 element of the sfermion mass matrix is — l/2rriQKsm 2# 2 3e ia . The magnitude of 
the FCNC between 2nd and 3rd generations is controlled by k sin 2^23, where #23 is the mixing 
angle in the unitary matrix. The phase parameter a also originates from the unitary matrix, 
and it will be the origin of a phase of the FCNC contribution. 

It is interesting that the flavor violation pattern in the lepton sector and the quark sector 
can depend on the SO(10) symmetry breaking vacua. Actually, in order to forbid a rapid proton 
decay, the quark flavor violation should be larger than the lepton flavor violation among the 
symmetry breaking vacua [23J. Namely, it is expected that kq, kuc, and Kpc are much larger 
than n L and he c - For example, if only the Higgs fields (8, 2, ±1/2) are light compared to the 
breaking scale (which is the most suitable case), one obtains kq = kjjc = k, D c, and only quark 
flavor violation is generated, while the lepton flavor violation is not generated. On the other 
hand, when the flavor violation is generated from the minimal-type of SU(5) vacua with type I 
seesaw, the quantities ft's have relations as k l ~ k D c, and kq, k uc , k E c ~ 0, effectively. Actually, 
when we take the threshold effect into account, it is expected that k,l is always larger than 
Kpc since the right-handed Majorana mass scale is less than the scale of colored Higgs mass. 
Therefore, the existence of b-s transition indicated by the experimental results in Fermilab 
predicts the sizable lepton flavor violation in the minimal-type of SU(5) model. Thus, if the 
results of large B s -B s phase is really an evidence of NP, the GUT models are restricted severely 
[H [9j [TDJ . Therefore, investigating the quark and lepton flavor violation is very important to 
decipher the GUT symmetry breaking especially when the B s -B s phase is large [9]. 

3 B s -B s mixing and direct dark matter detection 

Let us briefly see the phase of B s -B s mixing. We use the model-independent parameterization 
of the NP contribution: 

where 'full' means the SM plus NP contribution, M{f = Mf 2 M + Mf 2 p . The NP contribution 
can be parameterized by two real parameters Cb s and <pB a - The time dependent CP asymmetry 
(5* = sin0 s ) in B s — > J/ip(f) is dictated by the argument of Mflf : 4> s = — arg^f^ 11 ) an d thus (fi s = 
2(j3 s — (pB a )- It is important to note that large SUSY contribution is still allowed even though the 
mass difference of B s -B s [21] is kept fairly consistent with the SM prediction. This is because 
the mass difference, AMg s , can be just twice the absolute value of M\ \ . The consistency of 
the mass difference between the SM prediction and the experimental measurement just means 
Cb 3 ~ 1, and a large 4>b s is still allowed. For example, when Cb s — 1, the phase 4>b s is related 
as 2sin0 Bs ~ A^ P /A* M , where A^ P > SM = |M^ P,SM |. The argument of Mff, being free in GUT 
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models, is due to the phase in off-diagonal elements in SUSY breaking mass matrix (in the basis 
where Yd is a real diagonal matrix) , and one can choose an appropriate value for the new phase 
in the NP contribution. Therefore, the experimental data constrains A^ p /Af^, and therefore, 
K sin 2#23 is constrained for a given SUSY particle spectrum when the phase of B s -B s is large. 

In the MSSM with flavor universality, the chargino box diagram dominates the SUSY con- 
tribution to Mi2(B s ). In the general parameter space of the soft SUSY breaking terms, the 
gluino box diagram can dominate the SUSY contribution for a lower tan/? (i.e. tan/5 <^ 30). 
The gluino box contribution is enhanced if both left- and right-handed down-type squark mass 
matrices have off-diagonal elements [7], and therefore, it is expected that the SUSY contribution 
to the B s -B s mixing amplitude is large for the SO(10) model with type II seesaw, compared to 
the minimal-type of SU(5) model [9]. 

The box diagram does not depend on tan (3 (ratio of the vacuum expectation values of two 
Higgs fields) explicitly, whereas, the flavor changing Higgs interaction (through so-called Higgs 
penguin diagram) directly depend on tan /3, and the double Higgs penguin contribution to the 
B s -B s mixing can become more important than the box diagram when tan /3 is large and there 
is an off-diagonal element in the right-handed down- type squark mass matrix [T2j fl~3] . We 
note that the off-diagonal elements in the left-handed squark mass matrix is less important in 
order to generate a sizable double penguin contribution. This is because the chargino loop can 
generate the left-handed Higgs penguin contribution. Therefore, even in the minimal-type of 
SU(5) model, the double penguin contribution can be sizable when tan/3 is large. When the 
off-diagonal elements of left-handed squark mass matrix are generated, the left-handed flavor 
changing contribution to different processes (e.g., b — > 57) can be modified. 

The B s — > \m decay can be generated by a single Higgs penguin diagram [HI [13] . The decay 
amplitude is proportional to the muon Yukawa coupling, and thus the amplitude is proportional 
to tan 3 j3. Therefore, the branching ratio is proportional to tan 6 (3. Since it can be generated by 
a single left-handed penguin diagram, this decay occurs even in the universal SUSY breaking 
models like the mSUGRA (minimal supergravity) [25]. The current bound of the branching 
ratio is Br(5 s — > ////) < 4.7 x 10~ 8 [26]. When tan/3 is large, this bound gives an important 
constraint to the parameter space [HI EEl Ed E7J. In other words, one would expect that the 
B s — > decay will be observed very soon. 

When the lepton flavor violation is correlated to the flavor violation in the right-handed 
down-type squark as in the minimal-type of SU(5) model, the r — > fij decay will give us 
the most important constraint to obtain the large B s -B s phase [9j [10]. Furthermore, the 
squark masses are raised much more compared to the slepton masses due to the gaugino loop 
contribution since the gluino is heavier compared to the Bino and the Wino at low energy, and 
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thus the lepton flavor violation will be more sizable compared to the quark flavor violation. 
The current experimental bound of the branching ratio of r — > yy is [28] 




In order to allow for a large phase in the B s -B s mixing in the minimal-type of SU(5) model, a 
large flavor-universal scalar mass (often called mo) at the cutoff scale is preferable. The reasons 
are as follows. The gaugino loop effects are flavor invisible and they enhance the diagonal 
elements of the scalar mass matrices while keeping the off-diagonal elements unchanged. If 
the flavor universal scalar masses at the cutoff scale become larger, both Br(r — > yq) and 4>B a 
are suppressed. However, Br(r — > yy) is much more suppressed compared to 4>b s for a given 
k sin 2 $23 because the low energy slepton masses are sensitive to itlq while the squark masses 
are not so sensitive due to the gluino loop contribution to their masses. 

When tan (3 is large, the r — > yq constraint is relaxed for a large B s -B s phase, because the 
double-penguin contribution to the B s -B s mixing is proportional to tan 4 /? while the r — > yq 
is proportional to tan 2 /?. However, the B s — ► yy constraint becomes important in this case 
since it is proportional to tan 6 /?. As a result, the branching ratio of B s — > yy decay will have 
a lower bound in a given large tan (3 and SUSY spectrum when the phase of B s -B s mixing is 
large. This is because as follows: The double penguin contribution to the B s -B s mixing and 
the amplitude of B s — > yy are inversely proportional to m\, where is a CP odd Higgs 
mass. For a given tan/? and a large phase of B s -B s mixing, a larger k value is needed if rriA 
is supposed to be larger. Then, the parameter space is excluded by r — > yy constraint, due to 
the approximate relation Br(r — > yy) oc k 2 . Therefore, uia has an upper bound, and B s — > yy 
have a lower bound. 

Such constraints in the minimal-type of SU(5) have an impact on the neutralino dark matter 
which satisfies the recent WMAP result of relic density. The dark matter constraint is mostly 
satisfied in the minimal supergravity (mSUGRA) models by the following three scenarios. 

1. The stau- neutralino coannihilation region. 

2. The (nearly) Higgsino region (i.e. lightest neutralino has a large Higgsino component). 

3. The funnel region (i.e. the neutralinos annihilate through heavy Higgs bosons pole). 

At first, the stau-neutralino coannihilation region is not favored because the lighter stau is 
relatively light (almost degenerate with neutralino) in the region. When the stau is light, the 
r — > yy constraint will be severe and a large phase of B s -B s mixing is excluded. Secondly, the 
small Higgsino mass is not very favored either because it also enhances the r — ► yq amplitude 





(6) 
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unless the sleptons are very heavy. Such heavy sleptons do not explain the muon g — 2 anomaly 
[29J. The third solution (funnel region) is interesting in the current scheme to have a large 
phase of B s -B s mixing. In the funnel region, the lightest neutralino mass M^o is twice the 
mass of the heavy Higgs bosons (2M^o ~ tua)- In the models that we consider, the masses of 
heavier CP even Higgs boson and the CP odd Higgs boson are nearly degenerate. As we have 
mentioned, rriA has an upper bound for a given parameter space, and as result, the mass of the 
lightest neutralino is bounded from above. In the funnel region, m^ is almost fixed for a given 
gaugino mass. In the double penguin contribution of B s -B s amplitude and the B s — ► tua 
is a dominant parameter, and thus the branching ratio of B s — > fifj, can be predicted. 

The Higgs mass tua is also important for the spin-independent scattering cross-section since 
it gets the dominant contributions from the Higgs exchange diagrams. In much of the parameter 
space, the t-channel Higgs exchanges (h, H) dominate the proton-neutralino cross-section 0£o_ p . 
The spin- independent scattering cross-section can be written as [30J: 

a^o ~ -ml\{A u f u /m u + A c f c /m c + A*/ t /m t ) + (A d f d /m d + A'f./m, + A b f b /m b )\ 2 , (7) 

where, /, = (p\m q qq\p) /m p , and f u ~ 0.027; f d ~ 0.039; f s ~ 0.36; f c = f b = f t ~ 0.043 [31]. 
The down- and up-type quark Higgs amplitudes are 

A d,s,b = 92 m d,s,b f sina F h cos a F H \ 

AM W \ cos [3 m 2 h cos (3 m 2 H J ' U 



A u,c,t _ 9l m u,c,t ( cos a F h sin a F H 

AM W \sin (3 m 2 h sin (3 m 2 H J ' [) 

where g% is the SU{2) gauge coupling, the Higgs mixing angle a is usually small ~ 0.1 and 



F h = (N 12 - N u t&n6 w )(N u cosa + A^i 3 sina), (10) 
F H = (N 12 - N 11 t&n6 w )(N u sina - N 13 cosa). (11) 

Here Nu are the mixing amplitudes for the lightest neutralino Xi among the Bino (B), Wino 
(W) and the two Higgsinos (Hi, H 2 ): 

Xl = N n B + N 12 W + A^tfi + N 14 H 2 . (12) 

We see that the down-type quark Higgs amplitude, Eq.flS]), has larger contribution to the 
cross-section rather than the up-type ones A u,c,t due to larger f s . It is also apparent that the 
cross-section increases for smaller value of the heavy Higgs mass, m#, which scales with m^ 
and smaller values of fi which increases N13. 

Therefore, is also important for the proton-neutralino elastic scattering cross-section, 
and this is particularly the case for models with non- universal Higgs masses [32l [33| [34"] . Thus, 
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there is a correlation between the direct detection of the Milky Way dark matter and Br(B s — > 
/i/i) [17]. Furthermore, the dark matter cross-sections are also affected by its particle contents, 
i.e., whether the lightest neutralino is gaugino or Higgsino dominated. 

The current highest sensitivity of direct detection is about 5 x 1CT 8 pb for neutralino mass 
<^ 100 GeV [351 [36]. This is expected to increase to 2 x 10 -9 pb soon [37] for neutralino mass 
<S 100 GeV (for the neutralino mass we have considered in this paper the expected limit should 
be around 5 x 10~ 9 pb). For the models we consider, this constraint can exclude the parameter 
space with small rriA and/or small fi. 

4 Numerical results 

In order to illustrate the features described in the previous section, we plot the figures when 
the NP/SM ratio of the B s -B s amplitude is 0.5, A^ P /A$ M = 0.5, and the absolute value of the 
full amplitude is same as SM amplitude, Cb s = 1- Under these choices, one can obtain that 
|20 B J is about 0.5 (rad). We consider that the SUSY breaking Higgs squared masses, m 2 Hu and 
Tn 2 Hd , are not related to other scalar masses in order to make uia and [i free parameters, since 
these two parameters are important for Higgs penguin contribution and the proton-neutralino 
cross-section. 

Figure 1 is drawn in the case of minimal-type of SU(5) model. We choose the k values for 
the non-universality to be same for left-handed sleptons and right-handed down-type squarks, 
for simplicity. To draw, we choose tan/3 = 40, and the universal trilinear scalar coupling at 
GUT scale is zero (A Q = 0). The unified gaugino mass mi/ 2 is chosen to be 500 GeV and 800 
GeV, and the sfermion masses at GUT scale is chosen to be 500 GeV and 1 TeV. These ranges 
of mass parameters can be probed at the LHC. We plot By(B s — > /i/i) and the proton-neutralino 
spin independent cross-section (o" x _ p ) contours using black and green lines respectively. The 
blue strips are the 2-std region of the WMAP dark matter relic density, assuming that the 
entire dark matter content is made of neutralino. The gray shaded region is excluded by the 
experimental bound of Br(B s — > /i/i). The red shaded region (bottom-right corner) corresponds 
to the stau LSP and hence is disallowed by the dark matter requirement. The yellow shaded 
region is excluded by the experimental bound of Br(r — > /ry). As one can see from the figures, 
the t — > /i7 bound is relaxed for a larger sfermion mass mo- For a larger gaugino mass mi/2, 
the WMAP allowed strips for the funnel region (i.e., the vertical strips for 2M x o ~ m^) shift 
to the right. As mentioned in the previous section, the stau-neutralino coannihilation region 
(close to the red shaded region) and the small Higgsino mass (i.e. small /i) is excluded by the 
bound of Br(r — > /ry). On the other hand, the funnel region is still allowed. When mi/2 = 800 
GeV, the funnel regions (for fi < 1 TeV) shift to the region excluded by r — > fij. In order to 
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Figure 1: The — /i plane for minimal-type of SU(5) with non- universal Higgs masses, for 
tan/? = 40, A = and (m 1/2 ,m ) = (a) (500,500) GeV, (b) (500,1000) GeV, (c) (800,500) GeV 
and (d) (800,1000) GeV respectively. 
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tan P = 40, A = 0, m 1/2 = 800 GeV, m = 500 GeV 




1000 



1500 



1000 



Figure 2: The — /i plane for minimal type SO(10) models with non-universal Higgs masses, 
for tan (3 = 40, A = and (m 1/2) m ) = (a) (500,500) GeV, and (b) (800,500) GeV respectively. 

allow these regions, [i > 2 TeV and \i > 1.2 TeV are needed for m = 500 GeV and 800 GeV 
respectively. As a consequence, if /i is restricted to be less than 1 TeV, mi/ 2 is bounded and 
then the cross-section for the direct dark matter detection is bounded from below. Note that 
the neutralino mass, which depends mostly on mi/ 2 , in these plots, is large (about 200 GeV 
for (a) and (b), and 320 GeV for (c) and (d)). At these masses, the current sensitivity of the 
direct detection experiments (as can be seen in Fig. 4 of [35] and Fig. 4 of [36]) is still lower 
than what is needed to exclude more of the parameter space. With the expected increase of 
the sensitivity by an order of magnitude in the near future, this constraint would become more 
severe. 

Figure 2 is drawn for the case of SO (10) boundary conditions in the same way as for figure 
1. We choose the kappa values to be same for both left- and right-handed squarks. As is 
mentioned, depending on the SO (10) breaking vacua, the flavor non-universality in the slepton 
mass matrices can be reduced, and we choose them to be zero to escape from the r — > /i7 bound. 
With this choice, then, there is no upper bound for m^. Therefore, if Br(B s — > fifi) turns out to 
be small (< 10~ 8 ) and tan/5 large, we should adopt the SO(10) model. As we have mentioned, 
in the SU(5) case, when we take both the dark matter relic density and t —>■ fij bound into 
account, \i needs to be larger for a larger mi/ 2 , and consequently the direct detection cross- 
section would be too small to be probed at the ongoing XENON 100 experiment, as can be seen 
from figures l(c,d). In the SO(10) case, on the other hand, the direct detection cross-section 
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Figure 3: Correlation between Br(B s — > lili) and u x _ p for the minimal-type (a) SU(5) and (b) 
SO(10) cases, both for tan/3 = 40, A = 0, m 1/2 = 500 GeV and m = 500 GeV. 



even for larger values of (e.g. figure 2(b)) can still be large and could be probed very 
soon. This is, in part, because the stau coannihilation and the Higgsino dark matter solutions 
are also allowed along with the funnel region. If we find evidence for these solutions from the 
LHC, then the SO(10) model would be preferred. 

We comment that the boundary condition for the left-handed squark mass and right-handed 
up-type squark masses are different from the previous SU(5) boundary condition. Therefore, 
the RGE running for the SUSY breaking Higgs masses are different, and thus the stau LSP 
regions (red/pink shaded regions) for the SO(10) figures are different from the SU(5) ones. 

In Figure 3, we show the correlation between the Br(B s — > lili) and the proton-neutralino 
cross-section for SU(5) and SO(10) boundary conditions. The data points are picked up from 
the figure 1 and 2 for m = TO1/2 = 500 GeV, and are sliced for fixed values of the Higgsino mass 
li. The small circles represents the solution for the WMAP relic density. Since itla is almost 
determined for the funnel region of the WMAP solution and the li dependence of By(B s — ► lili) 
for a given phase of B s -B s mixing is not large, Br(5 s — > lili) is predictable as one can see from 
the figures. The proton-neutralino cross- sect ion, on the other hands, depends on the Higgsino 
mass hence varies with li. The current experimental sensitivity is still low, but if one of the 
solutions is correct, the direct detection of the dark matter should be possible in the near future. 
Note, however, that if the neutralino contributes only a part of the dark matter content (i.e., 
regions between two circles where the neutralino relic density is lower than the WMAP), the 
neutralino direct detection rate would be scaled down. In this case, we will need even higher 
sensitivity for the direct detection to exclude the parameter space regions. 
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5 Conclusion 



We have investigated the GUT models when the B s -B s mixing phase can become really large 
as indicated in the Fermilab experiments. We considered two cases: one is the minimal-type of 
SU(5) model with type I seesaw. The other is the minimal-type of SO(10) model with type II 
seesaw. The difference between the two boundary conditions is whether there exists a sizable 
off-diagonal element in the left-handed squark mass matrix. We emphasize that the sources of 
FCNC in the GUT models will be restricted if the large phase of B s -B s mixing persists in the 
upcoming result from Fermilab. 

In the case of a large tan/3, the double penguin contribution dominates the SUSY contri- 
bution to the B s -B s mixing when the pseudo Higgs mass is not too heavy. Especially in the 
minimal-type of SU(5) model, the pseudo Higgs mass should be low enough to satisfy the ex- 
perimental constraint from r — > /ry decay, and because of this, the branching ratio of B s — > 
is sizable and can be detected very soon. The r — > fi'y constraint also restricts the Higgsino 
mass \i and the slepton masses, and it may exclude some solutions of the relic density of the 
neutralino dark matter. In fact, the stau-neutralino coannihilation and the Higgsino dark mat- 
ter solutions are not favored for the large phase of B s -B s mixing in the minimal-type of SU(5) 
model. On the other hand, the funnel solution (in which the neutralinos annihilate through the 
heavy Higgs bosons pole) is favored in this case. For the funnel solution, the pseudo Higgs mass 
is almost determined, and the branching ratio of the B s —>■ decay is more predictive and 
can be observed soon for the values of soft masses which can be probed at the LHC. The direct 
detection cross-section depends on the heavy Higgs mass and the Higgsino mass, and correlated 
to the B s — > fifj, in the case of a large phase of the B s -B s mixing. The dark matter-nuclear 
cross-section could be in the range to be detected very soon in the upcoming experiments in 
both SU(5) and SO(10) models. 

In this paper, we have concentrated on the importance of the 2nd and 3rd generation FCNC 
effects such as Br(r — > fij) and <$b s correlation in GUT models, since they can be correlated 
directly by the 23 mixing. The constraints from Br(/i — ► ej) decay, K-K and B^-B^ mixings 
may also be important, but these effects depend on the details of the flavor structure which 
can have a freedom of cancellation. We refer to the Ref. [22j for an analysis of flavor violation 
including the first generation. 

Acknowledgments 

This work is supported in part by the DOE grant DE-FG02-95ER40917 (B.D. and Y.M.) and 
the NSF grant No. PHY-0652363 (Y.M.). 



14 



References 



[1] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 100, 161802 (2008) 
|arXiv:0712.2397l [hep-ex]]; V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 101, 
241801 (2008) |arXiv:0802.2"255l [hep-ex]]; D. Tonelli [CDF Collaboration], arXiv:0810.3229l 
[hep-ex] . 

[2] A. Lenz and U. Nierste, JHEP 0706 (2007) 072 |hep-ph/061216~7 ]; M. Bona et al. [UTfit 
Collaboration], larXiv: 0803 .0659 [hep-ph]; larXiv:0906.0953l [hep-ph]; L. Silvestrini, Nucl. 
Phys. Proc. Suppl. 185, 41 (2008); E. Barberio et al. [Heavy Flavor Averaging Group], 
larXiv:0808.1297l [hep-ex] . 

[3] F. Gabbiani and A. Masiero, Nucl. Phys. B 322, 235 (1989); J. Hagelin, S. Kelley and 
T. Tanaka, Nucl. Phys. B 415, 293 (1994); F. Gabbiani, E. Gabrielli, A. Masiero and 
L. Silvestrini, Nucl. Phys. B 477, 321 (1996) [hep-ph/9604387| . 

[4] F. Borzumati and A. Masiero, Phys. Rev. Lett. 57, 961 (1986); J. Hisano, T. Moroi, 
K. Tobe, M. Yamaguchi and T. Yanagida, Phys. Lett. B 357, 579 (1995) |hep-ph/950l"407| . 

[5] L. J. Hall, V. A. Kostelecky and S. Raby, Nucl. Phys. B 267, 415 (1986); R. Barbieri and 
L. J. Hall, Phys. Lett. B 338, 212 (1994) [hep-ph/9408~40~6| ; J. Hisano, T. Moroi, K. Tobe 
and M. Yamaguchi, Phys. Lett. B 391, 341 (1997) [hep-ph/ 96 05296| . 



[6] T. Moroi, Phys. Lett. B 493, 366 (2000) |hep^p h/0007328|. 

[7] B. Dutta and Y. Mimura, Phys. Rev. Lett. 97, 241802 (2006) |hep-ph/0607T47] . 

[8] J. K. Parry, Nucl. Phys. B 760, 38 (2007) [ hep-ph/ 05 10305] ; Mod. Phys. Lett. A 21, 2853 
(2006) |hep-ph/0608~l~92j ; J. K. Parry and H. h. Zhang, Nucl. Phys. B 802, 63 (2008) 
jarXiv:0710.5443l [hep-ph]]. 

[9] B. Dutta and Y. Mimura, Phys. Rev. D 78, 071702 (2008) [arXiv: 0805 .29881 [hep-ph]]. 

[10] J. Hisano and Y. Shimizu, Phys. Lett. B 669, 301 (2008) [ arXiv:0805.3327l [hep-ph]]; 
J. h. Park and M. Yamaguchi, Phys. Lett. B 670, 356 (2009) [ arXiv: 0809 .26141 [hep-ph]]. 

[11] B. Dutta and Y. Mimura, Phys. Lett. B 677, 164 (2009) |arXiv:0902.00T6l [hep-ph]]. 

[12] C. Hamzaoui, M. Pospelov and M. Toharia, Phys. Rev. D 59, 095005 (1999) 



|hep-ph/ 9807350|; M. Gorbahn, S. Jager, U. Nierste and S. Trine, larXiv:0901.2"065 [hep- 
ph]. 



15 



[13] A. J. Buras, P. H. Chankowski, J. Rosiek and L. Slawianowska, Nucl. Phys. B 619, 434 
(2001) [hep-ph/0107048| ; Phys. Lett. B 546, 96 (2002) |hep-ph/020724l| ; Nucl. Phys. B 
659, 3 (2003) |hep-ph/021014 5|. 

[14] S. R. Choudhury and N. Gaur, Phys. Lett. B 451, 86 (1999) |hep-ph/9810307] ; K. S. Babu 
and C. F. Kolda, Phys. Rev. Lett. 84, 228 (2000) [hep-ph/9909476| ; C. S. Huang, W. Liao, 
Q. S. Yan and S. H. Zhu, Phys. Rev. D 63, 114021 (2001) [Erratum-ibid. D 64, 059902 
(2001)] |hep-ph/0006~250 |. P. H. Chankowski and L. Slawianowska, Phys. Rev. D 63, 054012 
(2001) |hep-ph/0008~046 |; C. Bobeth, T. Ewerth, F. Kruger and J. Urban, Phys. Rev. D 
64, 074014 (2001) |hep-ph/0104284] ; Phys. Rev. D 66, 074021 (2002) |hep-ph/0204225] . 

[15] J. Foster, K. Okumura and L. Roszkowski, Phys. Lett. B 609, 102 (2005) 
|hep-ph/0410323| ; JHEP 0508, 094 (2005) |hep-ph/ 05 06146 ]; Phys. Lett. B 641, 452 
(2006) [hep-ph/0604l"2 1|. 

[16] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 180, 330 (2009) 
|arXiv:0803.0547l [astro-ph]]. 

[17] J. R. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, JHEP 0605, 063 (2006) 
[hep-ph/06031361 - 

[18] P. Minkowski, Phys. Lett. B 67, 421 (1977); T. Yanagida, in proc. of KEK workshop, 
eds. O. Sawada and S. Sugamoto (Tsukuba, 1979); M. Gell-Mann, P. Ramond and R. 
Slansky, in Supergravity, eds. P. van Nieuwenhuizen and D. Z. Freedman (North-Holland, 
Amsterdam, 1979); R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

[19] J. Schechter and J. W. F. VaUe, Phys. Rev. D 22, 2227 (1980); R. N. Mohapatra and 
G. Senjanovic, Phys. Rev. D 23, 165 (1981); G. Lazarides, Q. Shafi and C. Wetterich, 
Nucl. Phys. B 181, 287 (1981). 

[20] K. S. Babu and R. N. Mohapatra, Phys. Rev. Lett. 70, 2845 (1993) |hep-ph/9209215] . 

[21] B. Dutta, Y. Mimura and R. N. Mohapatra, Phys. Rev. D 69, 115014 (2004) 
|hep-ph/0402lT3| ; Phys. Lett. B 603, 35 (2004) |hep-ph/0406262 |; Phys. Rev. Lett. 94, 
091804 (2005) [hep-ph/0412105] ; Phys. Rev. D 72, 075009 (2005) lhep-ph/05073l"9] . 

[22] B. Dutta and Y. Mimura, Phys. Rev. D 75, 015006 (2007) |hep-ph/0611268| ; Phys. Rev. 
D 77, 051701 (2008) [ arXiv: 0708 .30801 [hep-ph]]. 

[23] B. Dutta, Y. Mimura and R. N. Mohapatra, Phys. Rev. Lett. 100, 181801 (2008) 
|arXiv:0712. 12061 [hep-ph]]. 



16 



[24] V. Abazov et al. [DO], Phys. Rev. Lett. 97, 021802 (2006) [hep-ex/ 0603029|; A. Abulencia 
et al. [CDF], Phys. Rev. Lett. 97, 062003 (2006) [hep-ex/0606027| . 

[25] A. Dedes, H. K. Dreiner and U. Nierste, Phys. Rev. Lett. 87, 251804 (2001) 
|hep-ph/0108037| ; R. L. Arnowitt, B. Dutta, T. Kamon and M. Tanaka, Phys. Lett. B 
538, 121 (2002) |hep-ph/0203069| ; J. K. Mizukoshi, X. Tata and Y. Wang, Phys. Rev. D 
66, 115003 (2002) [hep-ph/0208078) . 

[26] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 100, 101802 (2008) 
|arXiv:0712. 17081 [hep-ex]]. 

[27] J. R. Ellis, K. A. Olive and V. C. Spanos, Phys. Lett. B 624, 47 (2005) |hep-ph/ 05 04196] . 

[28] K. Hayasaka et al. [Belle Collaboration], Phys. Lett. B 666, 16 (2008) [arXiv: 0705 .06501 
[hep-ex]]; B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 95, 041802 (2005) 
[hep-ex/0502032] . 

[29] Muon g - 2 Collaboration, G. W. Bennett et al, Phys. Rev. Lett. 92, 161802 
(2004) [hep-ex/040T0~0~8| ; F. Jegerlehner and A. Nyffeler, Phys. Rept. 477, 1 (2009) 
|arXiv:0902.3360l [hep-ph]]. 

[30] G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267, 195 (1996) 
[hep-ph/9506380| . 

[31] J. R. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, Phys. Rev. D 71, 095007 (2005) 
|hep-ph/050200l] ; J. R. Ellis, K. A. Olive and C. Savage, Phys. Rev. D 77, 065026 (2008) 
|arXiv:0801.3656l [hep-ph]]. 

[32] E. Accomando, R. L. Arnowitt, B. Dutta and Y. Santoso, Nucl. Phys. B 585, 124 (2000) 
|hep-ph/0001019| . 

[33] J. R. Ellis, A. Ferstl, K. A. Olive and Y. Santoso, Phys. Rev. D 67, 123502 (2003) 
[hep-ph/0302032| . 

[34] H. Baer, A. Mustafayev, E. K. Park and X. Tata, JHEP 0805, 058 (2008) [arXiv:0802.3"384l 
[hep-ph]]. 

[35] J. Angle et al. [XENON Collaboration], Phys. Rev. Lett. 100, 021303 (2008) 
|arXiv:0706.0039l [astro-ph]]. 

[36] Z. Ahmed et al. [CDMS Collaboration], Phys. Rev. Lett. 102, 011301 (2009) 
|arXiv:0802.3"530l [astro-ph]]. 



17 



[37] E. Aprile, L. Baudis, for the XENON100 Collaboration, l arXiv: 0902 .42531 [astro-ph.IM]. 



18 



